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ABSTRACT 

Ultra high-energy cosmic rays (UHECRs) are beheved to be protons accelerated in magnetized 
plasma outflows of extra-Galactic sources. The acceleration of protons to ~ 10^" eV requires a source 
power L > 10**^ erg s~^. The absence of steady sources of sufficient power within the GZK horizon 
of 100 Mpc, implies that UHECR sources are transient. We show that UHECR "flares" should be 
accompanied by strong X-ray and 7-ray emission, and that X-ray and 7-ray surveys constrain flares 
which last less than a decade to satisfy at least one of the following conditions: (i) L > lO^^erg s~^; 
(ii) the power carried by accelerated electrons is lower by a factor > 10^ than the power carried by 
magnetic fields or by > 10'^ than the power in accelerated protons; or (Hi) the sources exist only at low 
redshifts, z <C 1. The implausibility of requirements (ii) and (Hi) argue in favor of transient sources 
with L > 10^°erg s^^ 

Subject headings: cosmic-rays — X-rays: general — gamma rays: observations — galaxies: nuclei 



1. INTRODUCTION 



The 



origin 



of ultra high-energy (> 10^^ eV) 
cosmic-rays (UHECRs) remains a myster y 
(jBhattachariee fc Sigl I [19981 : iNagano fc Watson I I2OOOI ). 
The sources have not been robustly identified, and the 
models of particle acceleration are challenged by the fact 
that the energy spectrum extends to > 10^° eV. Several 
observational clues suggest that the UHECR flux is 
dominated by extra-G alactic light nuclei: the sp ectrum 
flattens at ~ lO^^ eV (|Nagano &: Watson I [2000) . there 
is evidence for a composition change from he avy to 
light nuclei at ~ 10^^ eV (jBird et al. 1993; Abba si et all 
|2005[) . and the UHECR arrival direction distribution is 
nearl y isotropic (|Finlev fc Westerho3l2004l:rAbbasi et all 
I2004D . The recent detection of a weak anisotropy in 
the arrival distributio n of > 6 x 10^^ eV cosmic-rays 
(Auger Collaboration c t~an |2008[ ). is consistent with 
that predicted by assuming that the spatial distribution 
of UHECR sources co rrelates with the large-sca l e dis- 
tribution of galaxies ([Waxman. Fisher fc Piran ] [19971 : 
iKashti fc Waxmanl[2008li ! 

Although the identity of the UHECR particles is 
uncertain, we will assume here that they are pro- 
tons. This assumption is motivated by two arguments. 
First, the observed spectrum of cosmic-rays with en- 
ergies > 10^^ eV is consistent with a cosmological 
distribution of proton accelerators producing (intrinsi- 
cally) a power-law spectrum of high energy protons with 
d log N/d log E fti —2, for the number as a function 
of energy E (lWaxmanlll995bl : [Bahcall fc Waxmanl l200l 
IKashti fc Waxmanll2008n~ This intrinsic power-law spec- 
trum is consistent with that expected in models of par- 
ticle acceler ation in col lisionlcss shock s, for both non- 
relativistic (jBlandfold fc Eichlerl Il995f) and relativistic 
shocks (Waxman 2006; see however Keshet 2006). Sec- 
ond, the leading candidates for extra-Galactic sources, 
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namely gamma-ray bursts and active galactic nuclei (see 
below), are expected to accelerate primarily protons. 

Robust model-independent considerations imply that 
UHECR protons can only be prod uced by s ources with 
an exceedingly high power output (| WaxmanI ^004) . L > 
r'^P~^10'^^ erg s^^, where F and /3c are the Lorentz fac- 
tor and characteristic velocity associated with plasma 
motions within the source^. Since no steady source 
above this power threshold is known to exist within the 
100 Mpc GZK horizon, the distance to which the prop- 
agation of ~ 10^*^ eV protons is limited by the ir inter- 
action with the cosmic microwav e background (|Greisenl 
[1961 iZatsepin fc Kuz'minI [l969f) . the UHECR sources 
are most likely transient. A possible alternative is, of 
course, an unknown class of "dark sources", which pro- 
duce little radiation and therefore remain undetectable 
by telescopes. 

Only two types of sources are known to satisfy the 
above minimum power requirement: active galactic 
nuclei (AGN) - the brightest known steady sources, 
and gamma-ray bursts (GRBs) - the brightest known 
transient sources^. The absence of AGN with L > 
10^^ erg s~ ^ within the GZK horizon had motivated 



iFarrar fc Gruzinovl ()2008[ ) to suggest that UHECRs may 
be produced by a new, yet undetected, class of short du- 
ration AGN flares resulting from the tidal disruption of 
stars or accretion disk instabilities. 

We show in §[2]that if electrons are accelerated together 
with the protons in UHECR-producing flares, then their 
radiative losses will produce a bright flare of X-ray and 
7-ray photons. We then show in § [3] that existing X-ray 
and 7-ray surveys already put stringent constraints on 
the properties of UHECR flares. In § |4] we discuss the 
possibility of "hiding" the X-ray emission. Our conclu- 

* Somewhat more stringent li mits may be obtained b y specify- 
ing the acceleration process; see INorma.ii et al.l ([199D, IWaxmanI 
Ill995al) . and §0 below. 

^ It was recognized early on l IHillasI ([1984') and references therein) 
that while highly magnetized neutron stars may also satisfy the 
minimum power requirement, it is difficult to utilize the poten- 
tial drop in their electro-magnet ic winds for p roton acceleration to 
ultra-high energy (see, however. lAronsI II2003I )). 
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sions are summarized in § [5l 

Throughout our discussion, we consider a scenario in 
which the flare is associated with ejection of magnetized 
plasma from the source, and where the charged parti- 
cles are accelerated within the magnetized outflow. The 
non thermal emission from a wide range of sources is de- 
scribed within the framework of such a scenario. This 
includes AGN jets and GRBs, as well as t he transient 
AGN flares proposed by iFarrar fc Gruzinovl £2008). We 
parametrize the UHECR flares by their power, L, dura- 
tion. At, characteristic ejection speed /3c, and rate per 
unit volume in the local Universe, n. The fractions of 
the total energy output carried by protons, electrons and 
magnetic fields are denoted by e^, Ce and es, respectively, 
and we assume that the energy spectrum of accelerated 
electrons is similar to that of accelerated protons with 
a power-law index, d\ogN/d\ogE « —2. This index is 
expected for astrophysical sources which accelerate par- 
ticles in strong collisi onless shocks (jBlandfold fc Eichleil 
Il995t IWaxmaia l2006 fl. and is inferred from the radia- 
tion observed from a var iety of high energy sources, such 
as sup ernova remnant s (jRevnolds fc Elhson I [1993) and 
GRBs (|Waxmanlll99l . 

It should be pointed out that the composition of the 
jets of high energy sources is unknown. Two classes of 
models are generally discussed: jets where the energy 
flux is dominated by the plasma kinetic energy, and jets 
where most of the energy is carried by electromagnetic 
flux. For the "kinetic" jets, a plausible mechanism exists 
for energy dissipation, particle acceleration and radiation 
emission, namely internal collisionless shocks within the 
outflow. Within this mechanism, a particle distribution 
following d log N/ d log E' « —2 is naturally expected. For 
the "electromagnetic" jets, the mechanism of energy dis- 
sipation and particle acceleration is not well understood. 
We will assume that a d log A^/d log £' « —2 particle dis- 
tribution is generated in this model too, as suggested by 
observations. 

As explained at the end of § 12. 1[ our conclusions are 
valid for both spherical and jetted flows. Throughout 
the paper, L stands for the " isotropic-equivalent power" 
(i.e. for a flow which is conical rather than spherical, 
L stands for the power that would have been carried by 
the flow had it been spherically symmetric) , and h stands 
for the "isotropic-equivalent rate density" (i.e. the rate 
inferred under the assumption of spherically symmetric 
emission) . 

2. FLARE PROPERTIES 
2.1. Rates and luminosities 

We define a transient source to have an active phase of 
duration At shorter than the time delay AtcR between 
the photon and the UHECR arrival times. With this 
definition, a "steady" source is one which is still active 
when the UHECRs from it are being detected. The ar- 
rival time delay originates from deflections of the charged 
UHECRs by intergalactic magnetic fields, and can be ex- 
pressed in terms of the deflection angle, 9, and propaga- 
tion distance, d, as AtcR ~ O^d/Ac. The deflection angle 
is hmited to < l°(d/100 Mpc)^/^(-£;/10^° eV)"^ (se e de- 
tailed discussion in §2.2 of lKashti fc Waxmanll2008l ). and 
therefore 

Aicfl < 10''-5(d/100 Mpc)2(£;/102" eV)-2 yr. (1) 



For transient sources, the apparent number density 
of UHECR sources is energy dependent and given by 

flAtcR- 

The required number density of active flares is 
obtained from the observed energy production rate 
of UHECR protons per comoving volum e, e = 
E^dh^/dE = 0.7 ± 0.3 x IQ-^^erg Mpc-^ yr-i (jWaxmanl 
ll995bt[Bahcan fc WaxmiMl[200l . giving 

nAt = -4tt = 3.2 X lO-io^M^Mpc-3, (2) 
EpL/ A epLij 

where, £44 = e/10'''*erg Mpc~'^ yr"^, ^47 = 
(i/10^''erg s^^), epL is the total energy output in 
protons, and epL/K is the energy production per 
logarithmic proton energy {E) interval in the ob- 
served energy range. For the acceleration spectrum of 
strong collisionless shocks, E^drip/dE ^ const, we get 
A = ln(£;„iax/-Emin) and so Ai ee (A/10) - 1. 

The flare duration is limited by the absence of UHECR 
sources with multiple events, the so-called 'repeaters'. 
The absence of repeater s sets a lower limit on the number 
density of sources, n ([Waxman. Fisher fc Piran 1 fl997l : 
iKashti fc Waxm an 2008V This can be derived by not- 
ing that the nearly isotropic distribution of the ~ 30 
Auger events of energy > 6 x 10^^ eV requires tens of 
sources to be active within « 200 Mpc (the propaga- 
tion distance of protons with energies > 6 x 10^^ eV) 
within Auger^s field-of-view, implying n > 10~^'^Mpc~'^. 
For transient sources this requirement implies fiAt > 
10-5-^(Ai/Atci?)Mpc"^. For a GZK horizon distance 
of 200 Mpc corresponding to UHECR energy oi E ~ 
6 X lO^^eV, we get AtcR < 10^-^ yr, and therefore 

riAi > 3 X 10-i°7i_5 ( Mpc-3, (3) 
VlOyr/ 

where n_5 = (n/10"''^Mpc"^). Using Eq. (O we find 

Ai< lOnll^M^yr. (4) 

The number density of associated photon flares may 
be obtained by assuming that: (i) the accelerated elec- 
trons have the same initial power-law index for their 
energy spectrum as the protons, and (ii) the electrons 
lose all their energy to radiation (see W2.2\ below for the 
justiflcation of the latter assumption). The photon lu- 
minosity per logarithmic frequency (v) interval is then 
vLi, = eeL/2A — {ee/2ep){epL/ A), where the factor of 2 
is introduced since typically v <yi E^. This implies that 
the number density of active photon flares with a lumi- 
nosity > vLi, is 

nAi=:^^ = 1.6xlO-i"^-4^Mpc-3, (5) 
2ep vLn Ep (vLyj/^Q 

where {vLy)^^ = (i^Lj^/lO^'^erg s^^). 

Requiring that the acceleration time iacc be smaller 
than the plasma expansion time tdyn and the proton en- 
ergy loss time tin.ss , sets lower limits on L and the outflow 
Lorentz factor, F (|Waxmanlll995af ). In the following, we 
briefly describe these limits and derive the implied con- 
straints on the photon luminosity. Assuming that ac- 
celeration results from electromagnetic processes within 
the outflowing plasma, the acceleration time must exceed 
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the Larmor gyration time of the accelerated particle^, 
iacc > 2-KfRL/c = 2Tr fE'/eBc, where E' = E/T and / 
is a dimensionless factor of order a few which depends 
on the details of the acceleration mechanism, and where 
the various times are defined in the plasma rest frame. 
Requiring tacc < idyn = r/T/3c, where r is the radial 
distance from the source at which particle acceleration 
takes plac e, this implies B > fE//3er and equivalently, 
(jWaxman 19953) 

esL > 2 (^^) ' c = 6.6 X lO^'fjEi, erg (6) 

where ebL = Anr'^cT'^B'^ /8tt, and E20 = (E'/lO^o eV). 
The minimum photon luminosity is therefore 

iyL^>3.3 X lO^s^-^— £:|oerg s"^ 

>8.6 X iQ-^^L^E^oeig s-\ (7) 

One of our primary objectives is to demonstrate that 
exceptionally powerful flares with L > 10^° erg are 
required. In what follows we limit the discussion to 
flares with F < 10^'^, since a higher F implies L > 
lO^^e^^erg s^^. 

A lower limit on the bulk Lorentz factor F is 
set by requiring that the synchrotron loss time 
would exceed the acceleration time, iacc < ^loss = 
67rF(mpC^)^/ci3^f7T(me/mp)^i?, where irip and mg are 
the proton and el ectron masses. U sing r < 2F^cAt, this 
condition implies (|Waxmanlll995af ) 

\ one J \mp mpC^ J \ Zd" J 

(8) 

or numerically, 

F > l.lf^/^E%\esL,,Y'^^M-,^l\ (9) 

Hereafter, we drop the dependence on [3 since the flow 
is required to be at least mildly relativistic with /? « 1. 
Combining this result with Eq. ([G]), we get 

cbL > 7.9 X 10*V^-E2^^Aiy//2 gj.g g-i^ (;lo) 

The constraints on L and F are the same for a spher- 
ical and conical (jet-like) outflow as lo ng as the openin g 
angle of the jet 9j is larger than 1/F (|Waxmanlll995ar ). 
Thus, the constraints in equations (|B))- pi?l) apply in both 
cases, provided that L is interpreted as the isotropic- 
equivalent power. In the case of jets, there could be a 
discrepancy between the apparent number of UHECR 
sources and photon sources, in case the deflection an- 
gle of CRs by the intergalactic magnetic field is larger 
than max[0j,l/F]. For the Lorentz factors considered 
here, F < 10^'^, the magnetic deflections are smaller than 
1/F > 1/30 = 2° (see the opening paragraph of this sec- 
tion), implying that we should see the same sources in 
both photons and UHECRs. Under these circumstances, 
the results in equations ©-(O hold, provided that n 
and n refer to the isotropic equivalents quantities. This 
also implies that we can use isotropic equivalent lumi- 
nosities in the luminosity functions. 

® For acceleration in coUisionless shocks, iacc is larger than the 
Larmor time by a factor (c/f)'^ where v is the shock velocity in 
the plasma rest frame. 



2.2. High energy photon spectrum 

We next show that our estimate of the photon lumi- 
nosity, vLi, = eei/A, is valid for photon energies above 
~ 1 keV. We start our discussion by justifying the fast 
cooling assumption for the electrons. 

If electrons are accelerated through a process simi- 
lar to that of the protons, their acceleration time to 
a given energy would be similar to that of the pro- 
tons. The maximum electron energy would be lower than 
that of the protons, owing to their higher cooling rate. 
Equating the electron acceleration time to the radiative 
cooling time, ZTrfjerrieC^/eBc ~ rueC^ / crr^ecu, yields 
7g ^ eB /2Tif(jTU. Here 7e is the electron Lorentz factor 
in the plasma rest frame and u is the electromagnetic en- 
ergy density of the radiation together with the magnetic 
fields. The synchrotron photons emitted by the highest 
energy electrons carry an energy, TH{0.3jgeB /rriec) = 
0.3F(?ie^/27r/(TTTOec)i?^/u. Assuming that the electrons 
lose most of their energy to radiation, 

/li^syn.max « ^ McV. (11) 

/(ee + es) 

In the plasma rest frame, the electron Lorentz factor, 7c, 
at which the dynamical time is comparable to the syn- 
chrotron cooling time, is given by GTrmeC^ / o'tcJcB^ — 
r/Tc. Approximating r = 2F^C(5t, where St < At 
is the characteristic variability time within the flare, 
we get 7c = QttT^ Stm^c^ / aTesL. The observed en- 
ergy of the synchrotron photons emitted by electrons 
with a Lorentz factor 7c is hvc = 0.3hTj^eB/meC — 
5.4\/27r2efic"/2me5iFVa2(eBi)3/2, namely 

Z.«yn,c = 0.1F8(eBi47)-3/2^iy^ GHz. (12) 

Using Eq. ^ we obtain 

I^sy„,c < 0.1/-3p5^-3^^^^Qjj2_ (^3) 

Equations p3)) and pT|) imply that the electrons lose 
most of their energy to radiation and that synchrotron 
emission would lead to a flat spectrum, i^L^ = egL/A, 
from i^syn.max down to the larger frequency among I'syn.c 
and I'syn.min, the characteristic synchrotron frequency of 
photons emitted by the lowest energy to which electrons 
are accelerated. The value of /ii'syn.min depends on the 
details of the acceleration mechanism. As we show be- 
low, /it'syn.min «C 1 cV for acceleration in internal shocks 
within an expanding wind, and a flat spectrum, i^L^ — 
eg-L/A, is expected down to optical frequencies. However, 
the X-ray luminosity is expected to be comparable to 
eeL/ A also for /ii'syn.min Well abovc the X-ray band, ow- 
ing to electron cooling. Consider the extreme case where 
all electrons are accelerated to the maximum energy. 
In this case the luminosity is vL^, — eeL at ft.i'syn.inax, 
and the cooling electrons would produce a luminosity 

vLi,/eeL « (j^/j^syn,max)"'^^^ at ^syn.c < V < ^syn,max- 

This implies a hard X-ray luminosity, in the 10-lOOkeV 
band, of ^ eei/30, which is comparable to eeL/ A. 

By specifying the acceleration mechanism, we may de- 
rive an estimate for /ii'syn,!!!!!!. If the wind power is car- 
ried by the kinetic energy of the outflowing plasma and 
the plasma is heated through internal shocks (which also 
accelerate particles) within the outflow, then the char- 
acteristic temperature of the protons is ^ epTUpC^ . This 
follows from the fact that relative motions within the 



4 



Waxman & Loeb 



plasma rest frame are expected to be mildly, but not 
highly, relativistic. Consider, for example, two equal 
mass elements moving along the same directions with 
Lorentz factors Fi ^ r2 ^ 1. The Lorentz factor of 
th ese ma ss elements in their center of mass frame is 
•\/ri/r2/2, implying that a mildly relativistic relative 
motion is obtained unless their respective Lorentz fac- 
tors are very different. Thus, if electrons are coupled to 
the protons and carry a fraction ee of the energy density, 
then the lowest energy electrons will carry an energy of 
Ee,min ^ {f-e / (■p)mpC^ , giving 



and 



{)m{eJepf{eBL„Y'^T-^5t-^ eV, (14) 



< 3 X 10-" r^V^E^^^St^. (15) 



Using Eq. ([9]) we also have 



(16) 

The emission at a photon energy 3> 10 MeV is dom- 
inated by inverse Compton (IC) up-scattering of syn- 
chrotron photons. For coUisionless shock acceleration in 
internal shocks with Ee^min ^ (ee/ep)wpC^, the gamma- 
ray luminosity at photon energies 10 MeV is vL^ = 
min[l, ee/eB]eeL/A. If £'e,min > {^e/ ^p)mpC^ , the IC 
emission may be limited to photon energies 100 MeV, 
and may be shifted beyond the observable range (> 
0.1 TeV). For eg 3> es, the synchrotron luminosity is sup- 
pressed to vLy = ee{eB/^eY^'^L/ modifying the factor 
(ce/es) in Eq. ([7]) to (ce/es)^/^ and the factor (ce/cp) in 
Eq. dSD to {eeeBf/^/ep. 

The high energy, > 100 MeV, emission may be sup- 
pressed by pair production. A photon of high energy 
E.^ 3> meC^ may interact with lower energy photons, 
E'^ ~ T'^{meC^Y / E^, to produce e+e~ pairs. The opti- 
cal depth is T^^ — n^{E!^)a^jr/T, where n^{E'^) is the 
co-moving number density of photons with observed en- 
ergy E'^, and cr^-y is the e"^e~ annihilation cross-section. 
Using n^{E'^) w vL^ / A-kt^cTE'^ and the lower limit on 
F (Eq.El), we find 



'77 



< 10- 



46 



Ell'\eBU,f/" 



E-y 



(17) 



We therefore conclude that pair production may suppress 
the 100 MeV flux only for {vL^) > lO'^'^erg s^^. 

3. LUMINOSITY FUNCTION CONSTRAINTS 

Equation (O provides the number density of active 
flares required to account for the observed flux of UHE- 
CRs as a function of flare luminosity. Figure [1] com- 
pares this result with the cumulative number density of 
high luminosity AGN a t z < 0.2 in the e nergy bands 
of 0.5-2 keV (ROSAT, [Mivaii et al.l [200l . 17-60 keV 
fINT EGRAL. iSazonov et al.l 120071^. 15-195 keV (Swift 
BAT, 'Tueller et al."2008l), and > 100 MeV (EGRET, 
I Chian g & Mukhcrjec 1998). At the high luminosities un- 
der consideration, all the soft X-ray (0.5-2 keV) sources 
in the ROSAT survey are identified (including stars and 
X-ray clusters in addition to AGN), and all but one of 
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Fig. 1. — The number density of active flares, ncRjiare. = nAt, 
required to account for the observed flux of UHECRs (Eq. [Sj, com- 
pared to the cumulative number density of bright e xtra-Galactic 
sources at various energy bands: 0.5-2 keV (ROSAT . IMivaii et aEI 
l2000h . 17-60 keV (INTEGRAL. ISaSonov et al .11 20071 '). 15-195 keV 
(Swift BAT. ITueller et al.l 12008) . and > 100 MeV (EGRET, 
IChiang fc Mukherie3IT998l ) . The measured luminosity in different 
bands is converted to i/Li, assuming a photon index of -2 (consis- 
tent with the observed spectra). The solid segments of the curves 
represent the measured component of the local (z < 0.2) lumi- 
nosity function (LF), whereas the dashed segments of the curves 
represent the LF component which is inferred by measuring the 
number density of bright sources at higher redshift and then evolv- 
ing it to 2 ~ using the LF evolution with z as measured at lower 
i^Li/. The dotted segments of the curves represent the upper limit 
on the number density in the luminosity range where no sources 
have been observed . The GRB number density and luminosity (e.g. 
IGuetta et "al|[2005l) are shown for comparison. 



the Swift BAT sources are identified^. We also note 
that obscuration by a high column density of hydro- 
gen of soft X-ray sources can not dramatically alter the 
source number density since the AGNs selected in the 
hard or s oft X-ray bands have similar number densities 
(see also [Silverman et al.l I2007D . At the highest energy 
band, > 100 MeV, where the angular resolution is poor- 
est, source identification is incomplete, but the contribu- 
tion of unidentified sources can not lead to a significant 
change in the statistics of sources. In particular, EGRET 
had detected 60 high-latitude point sources that have 
not been identified, comp ared with the 44 hi gh-latitude 
sources identified as AGN ([Chiang et al.ll9 95^. The hard 
X-ray (17-60 keV, -195 keV) luminosity function (LF) 
shown in Fig. 1 is given by 



nhx = 10 



-11 



(18) 



It is important to emphasize that the luminosity func- 
tion constraints depicted in Fig. 1 refer to all the known 
bright sources on the sky, and so our conclusions do not 
apply exclusively to AGN flares, but to any other class 
of flaring sources. 

We first consider the case of near equipartition be- 
tween electrons and magnetic fields, Ee/es ^ 1- In this 
case, the flare luminosity should be i^L^, > lO'^^erg s~^ 

^ The identification of BAT sources is more complete than that 
of INTEGRAL sources because of the accurate positions provided 
by the Swift XRT follow-up. 
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(see Eq. [?])• It is obvious from Fig. [l] or from com- 
paring Eqs. ([5]) and pS]) . that for Ce/f-B ^ 1 the ob- 
served number density of sufficiently bright sources is 
much smaher than the density of active fLares required 
to account for the UHECR flux, unless Ce/f-p <^ 1. For 
pL^ > lO'^'^erg the number density of active flares is 
limited to < 10^^'^Mpc"'^, which implies based on Eq. (O 
that ep/ee > 10'^{vL,y)^j , and epL = K{vL^){ep/ ee) ^ 
lO^^Aierg s^^. For vL^, ~ 10*®erg s^^, the number den- 
sity of sources appears to be consistent with the required 
number density of active flares in Eq. (O for Se/ep ^ 0.1. 
However, the X-ray sources identified could be candi- 
date UHECR sources only if they are transient, and only 
a small frac t ion o f the sources observed are variable. 
iGrupe et al.l (pOOlh examined 113 bright ROSAT AGN 
on a time scale of ~ 6 yrs, and found that only 3 showed 
a factor of 10 variation over this time scale (all others 
varied by a factor less than 3). Similarly, Winter et al. 
(2008) compared XMM-Newton and Swift XRT observa- 
tions of 17 sources and found fractional variations of only 
a few tens of percent over ^ 100 days (see their Table 
12), suggesting that < 3% of all hard X-ray sources have 
a lifetime of < 10 years. This implies that the number 
density of X-ray sources variable on ^ 5 min (the typical 
integration time in the analysis of Grupe et al. 2001 is 
^ 300 s) to ~ 10 yr time scale is 

nhX.y^r - 3 X 10-13(j,i^)-g2.2 y[^^-\ (^g) 

Comparing with Eq. (O, this implies that UHECR flares 
must satisfy ep/eg > 500 and CpL > 10^°Aierg s~^ for 
(i^i,y)46 = 1. A similar constraint is obtained using 
EGRET'S LF. 

The requirement €pL > 10^°Aierg s~^ may be avoided 
if the magnetic field energy density is much higher 
than the electron energy density, eg/es <C 1. For 
ee/es < 10-2, the minimum flare luminosity is (Eq. [7]) 
I'Li, < lO^^erg s"^, and the required number density 
of active X-ray flares in Eq. ([5]) is consistent with the 
number density of variable X-ray sources in Eq. (fT9)) for 
Ce/ep ~ 1. The gamma-ray luminosity is suppressed 
by Ce/cB with uL^ < lO^^crg s^^, a range in which 
the number density of sources is poorly constrained by 
EGRET. Next, we consider the eg/es ~ 0.1 regime. 
Here the minimum flare luminosity is vL^, ~ 10^^ erg s""'^ 
and the required number density of active X-ray flares 
is consistent with the number density of variable X- 
ray sources for eg/cp < 10"^ and with EGRET's LF 
for Ee/ep < 10"^ (see Fig. [1]). As mentioned in § [21 
IC emission may be shifted above the observable range 
(> 0.1 TeV) in electromagnetically-dominated outflows. 
For Ee/eB ~ 0.1 and Cg/ep < 10^^ we get Cp/es > 10, 
which implies that the outflow can not be electromagnet- 
ically dominated. This, in turn, implies that the flares 
should be accompanied by observable gamma-ray emis- 
sion, and hence that eg/cp < 10^^ must be satisfied. 

The preceding discussion implies that for a fiare du- 
ration in the range Ihr < At < 10 yr, the requirement 
CpL > 10^°Aierg s~^ may be avoided only for eg/e^ < 
10~2 or eg/cp < 10~^. These constraints are likely to 
improve in the near future with new 7-ray data from the 
recently launched GL^ ST satellite^, and with proposed 

* http://glast.gsfc.nasa.gov/ 



X-ray telescopes such as EXISJ^ . Next, we consider the 
case of flares with Ai ^ 10 yr. Since there is little in- 
formation on source variability on such time scales, all 
observed sources are fiare candidates. For At > 100 yr, 
the active flare number density is > 10~^Mpc~^ (see 
Eq.[3]). At this density, the X-ray LF requires (see Fig.[T]) 
an X-ray flux vL^ < lO^^erg s~^, which implies through 
Eq. El that eg/es < 0.1. The EGRET LF requires ei- 
ther eg/cp < 10"'^ or that the IC gamma-ray emission be 
shifted outside the observable range, which may be pos- 
sible for flares that are electromagnetically dominated. 

There is one important caveat to the above con- 
straints. The required number density of sources is low, 
< O.lGpc"'', so that that no source is expected to be 
detected within a distance of ~ 1 Gpc, which is the 
GZK horizon of particles with E ~ 10^^ eV. This im- 
plies that snapshot surveys can not provide useful con- 
straints on the local {z ^ 0) number density of high 
luminosity flares. For this reason, the z ~ LFs shown 
in Fig. [1] are not measured directly at high luminosi- 
ties, i/L^ > lO^'^'^erg s~^. Rather, the number density 
of bright sources is measured at a higher redshift and 
the local number density is inferred from the evolution 
of the LF with z as measured at lower values of i/L^. For 
example, the number density of soft X-ray sources with 
vLjj > lO^^erg s~^ is measured to be '--^ 5 x 10~^^Mpc~'^ 
at z ~ 1 and inferred (but not measured) to be much 
lower than ^ 10~-^2Mpc~^ at 2; ~ 0, based on the LF 
ev olution measured at lower uL^, (see, e.g. Figure 5 
of lHasinger et al]|2005D . Long-term monitoring surveys 
offer much better prospects for constraining the source 
population than snapshot surveys. For example, if the 
flare duration is a few days, then a survey that lasts for 
a year can put constraints that are ~ 100 times better 
than a snapshot survey. Upcoming surveys, such as Pan 
STARRS^" are expected to provide relevant data soon, 
and planned surveys such as LSST^^ will provide better 
constraining power in the future. 

We can not exclude the possibility that the number 
density of flaring sources with lyL^, > lO^^erg s~^ does 
not decrease towards z ~ as fast as the number den- 
sity of lower vL^ sources, and remains at a level of 
~ 5 X 10~^^Mpc~'^, which is marginally consistent with 
that required for the local production rate per unit vol- 
ume of UHECRs. However, such a scenario is unnat- 
ural since it requires two coincidences: the flares must 
become a dominant source of energy output only at 
I'Li, > lO^'^^erg (or else they would modify the ob- 
served LF evolution at lower vL^,)^ and exist only at 
z ~ (or else we would observe them at z > 0.5). 

4. "DARK, PROTON-ONLY" FLARES 

It is difficult to rule out a scenario in which the 
UHECR flares involve " electromagnetically-dark" or 
"proton only" flares. Although there is currently no evi- 
dence or physical reasoning to motivate the consideration 
of a new class of hidden sources, we nevertheless discuss 
its required properties for the sake of generality. 

Since the cross-section for inelastic pp collisions is much 
smaller than the Thomson cross-section, ctt — 6.7 x 

^ http://exist.gsfc.nasa.gov/ 

http:/ /pan-Starrs. ifa.hawaii.edu/ 
http://www.lsst.org/ 
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10^^^ cm^, the X-ray emission may be suppressed (with- 
out affecting proton escape from the source) by postulat- 
ing the UHECR source to be embedded within an opaque 
plasma cloud of column density, ^ cr^^ ^ lO^'^cm^'^, 
which is optically thick to Compton scattering. If the 
outflow is jet-like and relativistic, scattering within the 
cloud would suppress the X-ray luminosity by a factor 
> r^/^j', where we assume that the jet opening angle 
9j > 1/r. A suppression of the expected X-ray luminos- 
ity, > lO'*^ erg s""'^, by a large factor, > F'*, would allow 
a sufficiently high number density of candidate UHECR 
sources to satisfy current limits on their electromagnetic 
luminosity. 

Similarly, since the cross-section for p7 (pion photo- 
production) collisions is much smaller than the Thom- 
son cross-section, the gamma-ray emission may be sup- 
pressed (without affecting proton escape from the source) 
by postulating the source to be embedded within an 
isotropic X-ray radiation field, with sufficiently high 
photon density to prevent the escape of gamma-rays 
through pair-production. The required photon column 
density, ~ lO^^cm^^, implies an X-ray luminosity of 
- 6 X 1043(^/xoi6cni)(£'^/lkeV)erg s-\ where R is the 
source size and is the energy of the background pho- 
tons. 

The predicted X-ray emission could also be suppressed 
by assuming that the source is embedded in an intense 
isotropic radiation field at IR, optical or UV frequen- 
cies, with an energy density far exceeding that of the 
magnetic field of the outflow, thus suppressing the syn- 
chrotron emission of the electrons by rapid IC cooling. 
For a relativistic outflow, the luminosity Liso associated 
with the isotropic radiation field could be much smaller 
than that associated with the outfiowing magnetic field, 
esL, since the energy density ratio in the plasma rest 
frame is F^Liso/esL. Note that such X-ray suppression 
does not change the conclusion of the second paragraph 
of § m that epL > lO^^erg s~^ is required for flares with 
X-ray luminosity of i^L^ > lO'^^erg s~^. The flux can be 
written as vL^, = (ee/ fx)L/A with a suppression factor 
/x > 1, implying that Eq. ^ should be modified to 



hAt 



2epfx vLi, 



1.6 X 10 



-10 



Mpc~^ 



(20) 

However, since the relation between proton and photon 
luminosities is also modified to epL = K(vLj^){fx^p/ ^e), 
the constraint on CpL is independent of fx- The X-ray 
suppression may affect, however, the constraint ee/es ^ 
1, that must be satisfied in the absence of X-ray sup- 
pression for flares with vL^ <C lO'^^erg s~^ (Eq. [7]). 
In the presence of X-ray suppression, we may write 
Eq. d?]) as F^ < 0.l{vLy)i^fx{^el<^B)~^- Combined with 
fx - T^Liso/esL < l(F/10)2Li,o,45 (see Eq.©, this im- 
phes Lisa > 10^^(ee/eB)(i^iiy)45^erg s~^. This isotropic 
luminosity requires the associated AGN to involve the 
most massive black holes in the Universe (~ lO^^M©) 
shining near their limiting (Eddington) luminosity, in or- 
der for the X-ray suppression to have a significant effect 
on our results. The absence of known sources of this 



extreme luminosity within the GZK horizon of UHE- 
CRs in the local Universe (see Fig. lib in Greene & 
Ho 2007, and Fig. 6 in Hopkins et al. 2007) rules out 
long-lived sources, but still allows for rare flares. We note 
that the minimum flaring time associated with the light 
crossing-time of the Schwarzscild radius of these black 
holes is > 1 day, and that the IR-UV variability of ac- 
tive AGN is observationally c onstrained to be weak on 
longer timesc ales jSc sar ct all 120061 120071 ) of up to sev- 
eral decades ()De Vries et alj|2005[ ). 

5. CONCLUSIONS 

The absence of steady sources of sufficient power to ac- 
celerate UHECRs within the GZK horizon of 100 Mpc, 
implies that UHECR sources are transient. We have 
shown that UHECR "flares" should be accompanied 
by strong X-ray and 7-ray emission. Figure 1 demon- 
strates that X-ray and 7-ray surveys constrain flares 
which last longer than ~ 5 min and less than a decade 
to satisfy at least one of the following conditions: (i) 
L > lO^^erg s~^; (ii) the power carried by accelerated 
electrons is lower by a factor > 10^ than the magnetic 
fleld power or by > 10^ than that carried by accelerated 
protons; or (iiij the sources exist only at low redshifts, 
z ^ 1 . The implausibility of requirements (ii) and ( Hi ) 
argue in favor of transient sources with L > 10^°erg s~^. 
The required luminosity is well above the brightest lu- 
minosity ever recorded in an AGN flare, and exceeds 
by two orders of magnitude the Eddington limit for a 
black hole of lO^^M©, the highest mass expected to 
exist within a distance of lOOMpc (|Lauer "etan[2007l : 
iNataraian fc Tr cister 2008). The results shown in Fig. 
1 exclude the regime of low flare luminosities considered 
by Farrar & Gruzinov (2008). 

The lower bound of 300s on the window of flare 
durations over which our constraints apply, originates 
from the integration time of the X-ray data used in Fig. 
1. For flare durations At < 300s, Eq. (jTU]) requires 
esL > 0.3 X 10^" erg s~^, not significantly different from 
the minimum luminosity inferred for longer flare dura- 
tions. 

We have also explored potential caveats to the above 
conclusions. Long-duration (> 100 years) flares which 
are electromagnetically dominated could evade the con- 
straints illustrated in Fig. 1 if their gamma-ray emission 
peaks outside the EGRET energy band. In addition, an 
unknown population of "electromagnetically-dark" flares 
is in principle possible (see 21 for details), although there 
is no physical motivation to make its existence natural. 

Future gamma-ray observations with GLAST and X- 
ray observations with EXIST would improve the statisti- 
cal constraints on the source population of UHECRs and 
potentially shed more light on their nature. 
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